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^ Abstract 

00 In a vortex-state magnetic nano-dislsPEl^ the static magnetization is curling in the plane, except in the 
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__i core region where it is pointing out-of-planeR^, either up or down leading to two possible stable states of 

^ opposite core polarity p. Dynamical reversal of p by large amplitude motion of the vortex core has re- 

cently been demonstrated experimentally^QSI^ raising fundamental interespHHI for potential application in 



^ magnetic storage devices^!. Here we demonstrate coherent control of p by single and double microwave 

53 pulse sequences, taking advantage of the resonant vortex dynamics in a perpendicular bias magnetic fielo^. 

S 

Optimization of the microwave pulse duration required to switch p allows to experimentally infer the char- 



acteristic decay time of the vortex core in the large oscillation regime. It is found to be more than twice 
shorter than in the small oscillation regime, raising the fundamental question of the non-linear behaviour of 
magnetic dissipation. 



Magnetic vortices are topological solitons with rich dynamical properties. The lowest energy 
excitation of the vortex ground state is the so-called gyrotropic modeP, corresponding to the gy- 
ration of the vortex core around its equilibrium position with a frequency in the sub-gigahertz 
rangeP'^. It is now established experimentally^ that the excitation of this gyrotropic motion 
leads to a dynamical distortion of the vortex core profile, as predicted by micromagnetic simu- 
lations and theoretical analysis. This distortion increases with the linear velocity of the vortex 
core and opposes the core polarity, until the critical velocity Vc ^ 1. 667^/^4^ (7 is the gyromag- 
netic ratio of the magnetic material and A^^ its exchange constant) is reached and the vortex core 
polarity is reversecP^. 

In zero magnetic field, dynamical control of the polarity is difficult due to the degeneracy of 
the gyrotropic frequencies associated to opposite polarities p = ±1, which can lead to multiple 
core switchin^^l^. Still, selective core polarity reversal is possible using a circularly polarized 
microwave magnetic field because the sense of the core rotation is linked by a right-hand rule to 
its polarit)^. Control of polarity switching can also be achieved by precise timing of non resonant 
magnetic field pulseP^, in a similar fashion as domain wall propagation in magnetic nanowire^. 

Resonant amplificatiorPH of the vortex gyrotropic motion enables to reverse the core polarity 
with minimum excitation powei*^^^^^, as it allows to concentrate the energy in a narrow frequency 
band. In this scheme, the damping ratio is an important parameter because it controls the minimum 
amplitude of the resonant excitation required to switch the core^"^. Here, it is shown that the 
damping ratio close to the reversal threshold is significantly larger than the one measured in the 
small oscillation regime. We associate this to the non-linear nature of the reversal proces^^^^^. 

Investigation of vortex core reversal using time-resolved imaging techniques based on X-ray 
magnetic circular dichroism is very powerful, as it enables direct determination of the vortex core 
trajector)ff^. However, it requires to average many events and the interpretation can become more 
complicated due to multiple reversal processes. In order to investigate the coupling between res- 
onant microwave pulses and vortex dynamics near the core reversal threshold, we use a magnetic 
resonance force microscope (MRFM) in combination with a bias magnetic field applied perpen- 
dicularly to the disk plane, that induces two distinct resonant gyrotropic frequencies associated to 
opposite core polaritie^^^^. This frequency splitting is used for a simple, single shot reading 
of the polarity state using MRFM, before and after microwave pulses have been applied to the 
vortex-state nano-disk. 

The MRFM detection setup^^, illustrated in Fig. la, is described in the Methods section. It 



allows to probe the magnetization dynamics in individual nanostructures"^, and in particular to 
spectroscopically determine the resonance frequency of the gyrotropic mode in a vortex-state 
nano-diskPI. The studied nano-disk with thickness 44 nm and diameter 1 /xm is made of NiMnSb 
alloy (see MFM image in Fig. lb), an ultra-low damping epitaxial material^^ (typical Gilbert con- 
stant a LLC is between 0.002 and 0.003). A gold antenna patterned on top allows to apply pulses 
of linearly polarized microwave magnetic field h in the plane of the disk. 

In our experiments, we choose the bias perpendicular field /igi^ = 65 mT in order to unambigu- 
ously discriminate the two different gyrotropic frequencies, or polarity states: resonant frequencies 
/_ = 217 MHz and /+ = 254 MHz respectively correspond to p = — 1 (core antiparallel to H, 
see Fig.lc) and p = +1 (core parallel to H, see Fig. Id). The microwave power employed to read 
the polarity state with MRFM is weak enough (P = — 19 dBm), so that p is not reversed during 
the reading sequencepl. We also emphasize that the bias perpendicular field ji^H = 65 mT is 
almost five times smaller than the static field required to switch the core polarit}^^^, so that no 
significant asymmetry between the two possible reversal processes (p=— ltop = +l and vice 
versa) is induced, as it will be clear from the data presented below. 

First, we explore in Fig. 2a the efficiency of single microwave pulses to switch the vortex core, 
depending on their duration, frequency and power. 11 _ pulses are defined as single pulses that 
reverse the polarity from p=— ltop = +l and 11+ pulses from p = +ltop=— 1 (see 
Fig.2b). The experimental data is acquired as follows: for n_ (11+) pulses, the p = — 1 (p = +1) 
state is first reset using an initialization pulse whose result is known to be fully deterministicpl. 
Then, a single microwave pulse of given duration w, frequency / and power P is applied, and 
the final polarity state is read using MRFM. An opaque coloured pixel marks pulse settings for 
which reversals are recorded with a 100% success rate, while a blank pixel means that no reversal 
is recorded. The contour plots presented in Fig.2a with different shades of red (blue) show the 
superposition of the results for n_ (11+) pulses with three different durations w: 100, 50 and 
20 ns. 

The frequency splitting introduced by the perpendicular bias field clearly appears in Fig. 2a. 
Due to the resonant character of the investigated switching process, a pronounced minimum in 
the power level required to reverse the vortex core is observed at a frequency /* (/+) close to 
the gyrotropic frequency associated to the initial core polarity, for each duration of the n_ (11+) 
pulses. This frequency discrimination introduced by H allows most of the recorded results to be 
fully deterministic. 



The minimum power P* required to reverse the core polarity increases as the pulse duration w 
decreases. We have plotted the dependences on w of the optimal frequency /^ (Fig. 2d) and of the 
optimal pulse energy £"* = P^'w (Fig.2e). In these plots, w ranges from 1 /xs down to 3 ns, with 
experimental points obtained from the analysis of data sets similar to those presented Fig.2a. There 
is a clear minimum in the optimal pulse energy that occurs around w = 50 ns. The position of 
this minimum yields an estimation of the characteristic decay time. If the pulse duration exceeds 
this characteristic time, the vortex core dynamics reaches steady-state before the end of the pulse, 
i.e., the pulse duration is uselessly too long. On the contrary, if the pulse duration is shorter than 
this characteristic time, some energy is wasted outside the resonance line. Due to the transient 
response of the vortex core to the excitation pulse, the optimal frequency /^ measured for pulses 
shorter than 50 ns shifts to higher frequenc)P^ (see Fig.2d), a general behaviour of any harmonic 
oscillator forced on time scales shorter than its decay time. 

To be more quantitative, we have performed a numerical calculation based on the analytical 
approach developed for zero applied magnetic fielcP^, with the new ingredient that due to the bias 
perpendicular field, /_ and /+ are distinct frequencies. The trajectory of the vortex core submitted 
to the microwave field pulse is calculated using Thiele's equation^'^, and Guslienko's criterion^ 
for core velocity is used to determine if the polarity is switched by the end of the pulse. The 
behaviour of the optimal pulse energy as a function of w can be reproduced with an analytical 
formula where the damping ratio rffoj-ced ^^ ^^ adjustable parameter (blue solid line in Fig.2e). 
The best fit is obtained for rffoj-ced = 0.018, corresponding to a characteristic decay time rforced = 
V(^forced2^/+) — 35 US. The predicted optimal pulse duration is Wmm = l-26rforced — 44 ns 
(see Methods section for details). Fig. 2c illustrates the good agreement of the model with the 
experimental data of Fig. 2a on the whole range of pulse parameters. 

The significant result here is that the damping ratio rffoj-ced = 0.018 inferred in the forced regime 
close to the core reversal threshold is more than twice larger than the value measured in the small 
oscillation regime of the gyrotropic mode, found to be rfvortex = 0.0075 (corresponding to a decay 
time Tvortex = 85 ns). The damping ratio dvortex is itself about three times larger than the Gilbert 
constant, measured from the linewidth in the perpendicularly saturated state of the NiMnSb disk 
and found to be aLLG = 0.0025 (see Supplementary Figure 1). The difference between rfvortex 
and aLLG corresponds to the expected increase due to topological renormalizatiorF^^. This en- 
hancement is produced by the strongly inhomogeneous spatial distribution of magnetization in the 
undistorted vortex state compared to the perpendicularly saturated state. Our experimental results 



of Fig. 2 give a quantitative estimation of the damping ratio rffoj-ced f^^ ^he vortex distorted by large 
oscillation^^^l^. In a further step, it would be interesting to evaluate how much of the observed 
increase of (ifoj-ced against rfyortex is due to topological renormalisation between the distorted and 
undistorted vortex states. As can be seen in Fig.2e, a practical consequence of this increase is to 
decrease the optimal pulse duration w^i^^ and to increase the corresponding minimum pulse energy. 

Next, we investigate in Fig. 3 the vortex core dynamics in the free decay regime by monitoring 
the core reversal induced by two consecutive microwave pulses. In these experiments, the pulse 
duration is set to t(; = 9 ns and the power to P = — 1.8 dBm, i.e., slightly below the minimum 
power P* = — 1 dBm required to reverse the core polarity p for such a short pulse. Contrary to a 
single pulse, two pulses separated by a delay r (Fig. 3a: n_-r-n_; Fig. 3b: U^-r-U^) can reverse 
p. The striking oscillatory dependence on the pulses carrier frequency and on the delay r observed 
in Figs. 3 a and 3b enlightens the phase coherent coupling between the vortex gyrotropic motion 
and the microwave excitation. During the pulse duration, the gyrotropic motion is forced at the 
pulse frequency, while in the free decay regime, the core oscillates towards its equilibrium position 
at its natural frequency, therefore acquiring a phase shift with respect to the excitation carrier. As 
a result, the efficiency of the second pulse to drive the vortex core to the reversal threshold, which 
is estimated out of ten attempts in Fig. 3, depends on the microwave frequency and on the delay 
between pulses in an oscillatory manner (the oscillation period scales as the inverse frequency 
detuning). 

To illustrate further this effect of coherence, similar experiments with a tt phase shift introduced 
between the two pulses have been carried out. It is clear that the regions where successful reversal 
are observed in Figs. 3c and 3d are complementary to those in Figs. 3a and 3b, respectively. Thus, 
phase control of the microwave excitation can trigger vortex core switching. 

Using the same approach as before, it is possible to calculate numerically such experimental 
phase diagrams. An excellent agreement with experiments is achieved, as can be seen in Fig.4a. 
In Fig.4b, we have plotted the calculated vortex core trajectory and velocity as a function of time 
corresponding to two consecutive 11+ pulses at / = 253 MHz separated by a 60 ns delay. In the 
top graphs, there is no phase shift between pulses and the time delay is such that the second pulse 
is efficient to amplify the gyrotropic motion from the beginning of the pulse, in contrary to the 
bottom graphs, where a tt phase shift is set between pulses. As a result, the vortex core is not 
reversed by the end of the second pulse in the latter case, while it is in the former. The obtained 
data sets also allow us to fit the free decay time, and the best agreement with the experiment is 



obtained for rfree — 53 ns. Analysis of the difference of this value with rforced should further shed 
some light on the precise non-linear nature of magnetic dissipation close to the reversal threshold 
and on its dependence on the amplitude of the vortex core motion. 
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METHODS 

A. Sample preparation. 

The magnetic nano-disk (thickness 44 nm, diameter 1 /xm) was patterned by standard e-beam 
lithography and ion-milling techniques from an extended film of NiMnSb grown by molecular- 
beam epitaxy on an InP(OOl) substrate^^. A 50 nm thick Si3N4 cap layer was deposited on top 
of the disk for protection and a broadband coplanar microwave antenna (300 nm thick Au) was 
subsequently evaporated on top of the patterned disk. The width of the antenna constriction above 
the disk is 5 /xm (see Supplementary Figure 2). 



B. Microwave setup. 

Injecting a microwave current from a synthesizer inside the antenna produces an in-plane lin- 
early polarized microwave magnetic field /i, oriented perpendicular to the stripe direction. In order 
to apply two consecutive pulses separated by a variable delay r , the carrier of the cw excitation is 
split into two branches, each of them being gated by independent mixers before being recombined. 
A phase shift can be introduced between the two branches using a delay line. The typical rise and 
fall times of pulses are 0.8 ns and the control on the delay r is better than 0.1 ns. The calibration 
of the amplitude of h yields the value /xq/i = 1.05 mT with an error bar of 25% for a dBm input 
power in the antenna. Therefore, the microwave power range in Fig. 2a corresponds to microwave 
field amplitude ranging from 0.2 to 0.7 mT, the power in double pulses experiments of Fig.3a to 
fioh = 0.85 mT, and the power employed to read the polarity to fioh = 0.1 mT. 



C. Initialization pulse. 

A n+ pulse with settings (w = 50 ns, / = /+, P = —11 dBmj^is used to set the initial 
polarity state to p = — 1. In fact, such an initialization pulse does not affect the p = —1 state 
and it transforms p = +1 into p = —1 (see Fig. 2a). Similarly, a n_ pulse (w = 50 ns, f = f-, 
P = —11 dBm) is used to set the initial polarity state io p = +1. 

D. MRFM detection. 

The MRFM setuppl is located inside a vacuum chamber (10"^ mbar) and operates at room 
temperature. The cantilever is an Olympus Biolever (spring constant k o^ 5 mN/m) with a 800 nm 
diameter sphere of soft amorphous Fe (with 3% Si) glued to its apex. MRFM spectroscopy is 
achieved by placing the centre of the magnetic spherical probe above the centre of the NiMnSb 
nano-disk. The separation between the sample and the probe is 5 = 1.5 /xm. The probe senses the 
dipolar force F^ proportional to the perpendicular component M^ of the magnetization of the nano- 
disk (see Fig. la). Ferromagnetic resonance spectra are obtained as a function of the microwave 
excitation frequency at a fixed bias field H. The microwave modulation is a cyclic absorption 
sequence, where the microwave power is switched on and off at the cantilever resonance frequency, 
fc c^ 11.8 kHz. The MRFM signal originates from the cyclic diminution of M^ of the nano-disk 
synchronous with the absorption of the microwave fielcP^^. The resulting force modulated at the 
mechanical resonance of the cantilever force produces a cantilever vibration amplitude enhanced 
by its quality factor Q c^ 4000, that is optically detected. 

E. Numerical calculations. 

The calculations presented in Figs.2 and 4 are based on linearised Thiele's equation^^, an ef- 
fective equation of motion for the vortex core position X in the disk plane: 

-GxX - DX + kX + fi[z xh]=0. (1) 

The first term is the gyroforce (G is the gyro vector), the second one is the damping (d = —D/\G\), 
the third one is the restoring force (k. is the stiffness coefficient), and the last one is the Zeeman 
energy with the spatially uniform external field h (z is the unit vector normal to the disk plane 
and details on fi and previously mentioned coefficients can be found in refs.'^ES). The resonant 



frequency of vortex core gyration is uug = i^/\G\. From Eq|T]one can find the instantaneous po- 
sition X{t) and the velocity v{t) of the vortex core in the presence or absence of the harmonic 
excitation field /i(t), as derived in the supplementary documents of refM If Guslienko's criterion 
on critical velocity^ is met at some time of the simulated pulse sequence (v{t) > K). the vortex 
core is assumed to have switched. Experimental sample size, gyrotropic frequencies /_ and /+ as- 
sociated to each polarity, and non-linear redshift of the frequency (3%) are used in our calculation. 
The minimal excitation amplitude hc{uj,w) to reach Vc by the end of the pulse depend^ on the 
excitation frequency uj and on the pulse duration w. Minimization of the energy E oc whl{uj,w) 
with respect to uu yields the dependence on w of the optimal frequency /^ plotted in Fig. 2d and of 
the optimal pulse energy E* plotted in Fig.2e: 

E*iw) = w'i^ = f ^y :^ ^. (2) 

In this expression, a is an experimental conversion factor between the input power in the an- 
tenna and /i^ determined from calibration. Therefore, two independent parameters are used to 
adjust the calculation with the data in Fig.2e: the damping ratio d to fit the overall shape (position 
of minimum) and the the critical velocity Vc to fit the absolute value of the energy. The opti- 
mal pulse duration Wj^i^ corresponding to the minimum energy E^^^ follows from Eqj2| w^nin = 
1.26/((ifoj-^g^cjG) = 1.26rforced (for the blue sohd line in Fig.2e, rf^^.^^^ = 0.018, rforced = 35 ns, 
hence E^^ is reached for w^i^^ = 44 ns). The fitted critical velocity in our NiMnSb disk is 
Vc ^ 190 m/s, in good agreement with the expected value^"^ of 225 m/s (the exchange constant of 
NiMnSb is A^^ = 6 pJ/mP^. The same parameters are used in calculations of Figs. 2 and 4. In 
Fig.4, the only new fitting parameter is the characteristic time rfree of the vortex core free decay 
(the relaxation is assumed to be exponential). We also allow a fine adjustment of the phase shift 
between the two pulses in order to get the best agreement with the data (25° in the "zero phase 
shift" experiment and 190° in the "tt phase shift" experiment; these small differences are ascribed 
to an imperfect delay line in the experimental pulse setup). 
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FIG. 1. Experimental setup and sample, a, A magnetic resonance force microscope (see Methods) is used 
to probe the vortex core dynamics of an individual vortex-state NiMnSb disk (diameter 1 fim, thickness 
44 nm). A soft cantilever, with a spherical magnetic probe attached at its end and placed at a distance 
5 = 1.5 fim from the sample, detects mechanically the vortex core dynamics. The bias magnetic field H is 
applied perpendicularly to the disk plane whereas microwave pulses of linearly polarized field h are applied 
in the plane, b. Zero-field magnetic force image of the magnetic vortex, where the dark spot at the disk 
centre reveals the core prepared in the p = +1 (c) polarity state. The opposite configuration is the ^ = —1 
(d) polarity state. 

^^ Ritchie, L. et al Magnetic, structural, and transport properties of the heusler alloys Co2MnSi and 
NiMnSb. Phys. Rev. B 68, 104430 (2003). 
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FIG. 2. Vortex core resonant reversal by a single microwave pulse, a, Efficiency of the n_ and n+ 

pulses (see panel b) to reverse the vortex core at jj^qH = 65 mT measured as a function of power and 
frequency (stepped by increments of 0.3 dBm and 1.2 MHz). For each pulse type, a coloured pixel (n_: 
red shade, 11+ : blue shade) marks a successful reversal. The transparency gives the switching probability 
averaged over 16 attempts. Experiments corresponding to three values of the pulse duration w are displayed. 
For each w we define (/^,P*), the optimal working point of the 11+ pulse located at the bottom of the 
corresponding contour plot (see ■). c. Numerical calculation (see Methods) of the experiments presented 
in a. This calculation is not valid within the shaded area, where multiple vortex core reversals can occur. 
d, e. Experimental (■) and calculated (lines) dependences on w of the optimal frequency /+ (d) and of the 
optimal pulse energy £^* = P*it; (e). The experimental points are obtained from the analysis of data sets 
similar to those presented in a, where w is varied from 1 fis down to 3 ns (the three ■ close to the minimum 
energy in e are inferred from those displayed in a ). The best agreement is obtained for a damping ratio 
^forced ^ 0.018, i.e., a characteristic decay time Tforced = 35 ± 4 ns. The absolute value of the energy is also 
fitted in the calculation, which enables to extract the critical velocity for vortex core reversal, Vc — 190 m/s 
(see Methods). Error bars on /+ and E* are absolute minima and maxima resulting from experimental 
uncertainties in a on the optimal working point ( /T , P*) associated to each pulse duration. 
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FIG. 3. Oscillatory dependence on frequency and delay between two consecutive pulses of the vortex 
core reversal efficiency, a-d, Number of switching events out of ten attempts as a function of the delay r 
separating the two pulses (it; = 9 ns, P = — 1.8 dBm) and of the carrier microwave frequency (stepped by 
increments of 3 ns and 4 MHz, respectively). The bias magnetic field is fioH = 65 mT. The initial polarity 
state is p = —1 in the left graphs (a, c) and p = +1 in the right graphs (b, d). As depicted in the right-side 
panels, the phase difference between the two pulses is zero in the upper graphs (a, b) and n in the lower 
graphs (c, d). 
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FIG. 4. Phase coherent control of vortex core reversal, a, Numerical calculation (see Methods) of the 
double pulse sequences presented in Fig. 3. The best agreement is obtained for a characteristic decay time 
Tfree = 53 ± 6 US in the free regime, b, Associated vortex core trajectory (left) and velocity (right) vs. time 
plotted for two 11+ pulses with settings r = 60 ns and / = 253 MHz (see pixel in a). For these settings, 
the vortex core is reversed when the phase difference between the pulses is equal to zero (top, see star) and 
not reversed when it is equal to tt (bottom). 
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FIG. 5. (Supplementary Figure 1) Damping ratio in the small oscillation regime, a, Frequency- 
perpendicular field dispersion relation of the uniform mode in the saturated state (fJ^oH > 0.8 T) and of 
the gyrotropic mode in the vortex state (^o^ < 0.8 T). Dashed lines are analytical expressions (Eqs. 
1 and 2 in ref.l^. b, Dependence on frequency of the damping ratio obtained from MRFM spectro- 
scopic measurements in the small oscillation regime (the typical power level is —30 dBm, corresponding to 
fj^oh = 0.03 mT). Dashed lines are guides to the eye. Lorentzian fits of the resonance peaks (inset shows 
such a fit at a fixed field in the vortex state) yield the centre frequency / and the frequency linewidth A/ 
together with the associated error bars, hence the damping ratio d = A// (2/). The damping ratio in the 
saturated state is equal to the Gilbert constant d = aLLG = 0.0025. In the vortex state, the damping ratio is 
renormalized by topology (Eq.l3b in ref.l^ZI). Experimentally, it increases from about 0.005 to 0.008 as the 
frequency decreases from 450 MHz down to 190 MHz (i.e., the perpendicular field decreases from 0.65 T 
down to —0.15 T). Close to i7 = 0, the experimental value d = dvortex = 0.0075 ± 0.001 corresponds 
to a decay time Tyortex = l/(^^/) = 85 ± 12 ns and is in good agreement with the predicted valuePl 
dyortQx = o^LLG [1 + l^ (^/^)/2] — 0.007 (R = 500 nm and 6 2^ 15 nm are respectively the disk and vortex 
core radii). 
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FIG. 6. (Supplementary Figure 2) Microwave antenna and NiMnSb disk sample. Scanning electron 
microscopy (SEM) image of the microwave antenna where the red rectangle points out the 5 fim wide 
constriction under which the NiMnSb disk (diameter 1 fim, thickness 44 nm) is located (see insert). The 
microwave antenna is wire bounded to a coplanar waveguide on the right side. 
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